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Abstract 

I construct a new universal differential equation (|^, in the sense of Rubel. That 
is, its solutions approximate to arbitrary accuracy any continuous function on any 
interval of the real line. 


1 Introduction 

[Rubel 1981] proved the surprising theorem: 

Theorem 1 

Given any continuous function ^ : M — M and any positive continuous function e : M — 
M_,_, there exists a C°° solution y of 

3y y y -4y y y +6y yyy + 

24yVV"-i2yVy"^-29yV'V"^+i2/" = 0 (r) 

such that \y(t) — (l)(t)\ < £(t) Vf G M. 

Note that Rubel’s function y is C°° but not real-analytic, typically having a countable 
number of essential singularities. The solutions of (R) have the additional properties that 
they may be taken to be monotone if (j) is monotone, and can be made to agree with (j) at 
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a countable number of distinct abscissae. We may rephrase the result by saying that all 
continuous functions are uniform limits of sequences of solutions of (R), or that solutions 
of (R) are dense in the space C(M) of eontinuous funetions. 

By a universal ODE I will mean any nth order homogeneous non-trivial polynomial 
ODE P{y,y,y ',... = 0 (that is, an algebraic differential equation or ADE) whose 

solutions have the same property as the function y in Rubel’s theorem. Universal differ¬ 
ential equations are of interest in the theory of eomputability by analog eomputers. Eor 
discussions of this field, see [Boshernitzan & Rubel 1985, Rubel 1983, Rubel 1992, Rubel 
1996]. 

There have been several recent developments, which I now describe. [Duffin 1981] 
found analogous universal differential equations whose solutions are C" for finite n> 3: 
Theorem 2 

The differential equations 

n^y"'y'^+3n{l—n)y'y''y+{2n^ — 3n + l)y'^ = 0 (Dl) 

and 

ny'^'y^ + {2 — 3n)y"'y'y'+2{n — l)y'^ = 0 (D2) 

are universal. 

(Note that [Duffin 1981] eontains several errors whieh I have corrected in the above state¬ 
ment of (pT]).) 

[Boshernitzan 1986] gave a general survey of the field and proved a new result for 
finite intervals without a smoothness assumption: 

Theorem 3 

There exists an ADE (of order less than 20), the polynomial solutions of which are dense 
in C{I) for any Unite interval I. 

The ADE was not explicitly constructed. Of eourse, by Weierstrass’ classieal theorem 
on approximation by polynomials, we know that the set of all polynomials P{I) on an 
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interval I is dense in C(/); the interest in theorems such as the above is that a subset of 
P{I) is still dense. 

[Eisner 1992] proved a new result for Lipschitz functions, and [Eisner 1999] con¬ 
structed a universal functional equation. (This paper also contains an error: the statement 
of Theorem 1.1 should say that the greatest coefficient is 39813120000, not 35831890000.) 

2 The new universal differential equation 

Theorem 4 


//// /2 T ////// I 0/1 -2\ fl'i A /r>\ 

y y -3y y y+2{l-n )y =0 (B) 

is universal for n > 3. 

PROOF: I will make use of the Jacobian elliptic functions which depend on a real param¬ 
eter me [0,1) and satisfy 


sn'(x, 

m) 

= 

cn(x, m) dn(x, m) 

cn!{x, 

m) 

= 

— sn(x, m) dn(x,m) 

dn'(x, 

m) 

= 

—msn(x, m) cn(x,m 

sn^(x, 

m) 

+ 

cn^(x,m) = 1 

dn^(x, 

m) 

+ 

msn^(x, m) = 1. 


These properties alone will suffice for the proof. I let g(ji;) = cn"(a:, m) and consider the 
expression 

ng"' (x) g^ (x) + bg” {x)g' {x)g{x) +cg'{xf. (*) 

We wish to find (if they exist) values of m, n, h, c which make the expression (^ identically 
zero. After substitution of the expression for g(x), (^ simplifies to 

(1—m) \2—b+{b+c — l>)n+n^'\ + 

(2m— \)n(b+c-\-n) ex?{x,m) — 
m \2-\-b-\-{b-\-c-\-'i)n-\-n^'\ cn"^(x,m) 
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and we have the solutions 


m = 0, b = 2 — 3n, c = 2(n —1); 
m=l/2, b=—3n, c = 2{n — \/n); 

m=l, b~—2 — 3n, c = 2(n + l). 


The first solution in faet reduees to (PID; the last is exeluded as the elliptie funetions are 
not defined for m = 1, leaving only the ease m = 1/2 to eonsider. The funetion g and its 
first n — 2 derivatives vanish when x = ±K{m), where K is the eomplete elliptie integral 
of the first kind. Thus, if we extend g to the whole real line by making it zero outside 
[—K{m),K{m)], the extended funetion is still a solution of (^. Also, affine transforma¬ 
tions of g of the form g{ax+l3) are also solutions by homogeneity, and so are sums of 
sueh funetions on disjoint intervals. Thus 



forms an S-module in the sense of [Rubel 1981]on [K{m),K{m)]. (That is,it is monotonie 
and its first n—l derivatives vanish at the endpoints.) Our eomplete solution y is then 
obtained by pasting together suffieiently many sueh S-modules. ■ 

Note that (pT|), (P2[), and (0) all reduee to the same ADE in the limit n —> oo. It would 
be interesting to know whieh other ADEs of the general form of (^) are universal. 
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